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\0 ' Abstract 

O '. 

Q>^ , We have applied the Meson Cloud Model to calculate inclusive momentum spec- 

^^ . tra of pions and kaons produced in high energy proton-proton and proton-nucleus 

o ■ collisions. For the first time these data are used to constrain the cloud cut-off pa- 

Q_i' rameters. We show that it is possible to obtain a reasonable description of data, 

(— I I especially the large xp {xp > 0.2) part of the spectrum and at the same time de- 

scribe (partially) the E866 data on d — u and d/u. We also discuss the relative 
^ . strength of the irN and vrA vertices. We find out that the corresponding cut-off 

j^ ' parameters should be both soft and should not differ by more than 200 AleV from 

each other. An additional source (other than the meson cloud) of sea antiquark 
asymmetry, seems to be necessary to completely explain the data. A first extension 
of the MCM to proton nucleus collisions is discussed. 
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1 Introduction 

Recent measurements |l], § have established the asymmetry in the distributions of up 
and down quarks in the nucleon sea, a result which can not be understood in terms 
of perturbative QCD. The presence of pions in the nucleon can naturally account for 
the excess of d over u. The role of mesons in deep inelastic scattering (DIS) was first 
investigated by Sullivan [Q. He suggested that some fraction of the nucleon's anti-quark 
sea distribution may be associated with the pion cloud of the nucleon. Several works 
developed this idea and gave origin to the Meson Cloud Model (MCM) [^ |^, ^. They 
are all based on the notion that the physical proton (p) may be expanded in a sum 
of virtual meson-baryon {MB) states. The probabilities of these states are not known 
apriori. They are commonly related to the probability of the splitting p —>■ MB, which, 
in turn, is calculated with a simple Feynman diagramm of meson emission. In these 
calculations the assumption is made that the proton is already an extended object and 
a form factor is assigned to the meson emission vertex. This form factor contains a cut- 
off parameter which must be adjusted by fitting experimental data. In some works it 
is adjusted to correctly reproduce the Gottfried Sum Rule (GSR) violation. The most 
consistent procedure, however, is to fix the cut-off by simultaneously analyzing data on 
hadronic collisions and parton distribution functions. This is the atitude adopted, for 
example, in 0. The cut-off choices vary over a wide range according to the experimental 
source used in the determination procedure. In ref. § a detailed discussion of the subject 
is presented. 

In a very recent analysis of the E866 data [Q], the authors used the MCM to fit the 
d{x) — u{x) as a function of x. The dominant intermediate MB states are vrA^ and vrA. 
The conclusion was that the cut-off associated with both states must be soft and also 
AvrAf > AttA- In particular, with dipole form factor A^^j^ = 1-0 GeV and A^a = 0.8 
GeV. A similar conclusion was reached in |Q]. With such choices it was possible to 
find "nearly exact accord" with experiment. This finding on one hand increases the 
confidence in the virtual meson-baryon picture but on the other hand imposes constraints 
on the MCM description of other high energy collision data. In many previous studies 
AjrAf and A^a were assumed to be roughly equal 0, |^. Shortly after the E866 analysis, 
Melnitchouk, Speth and Thomas ||T0| presented a new calculation of the d — u asymmetry 
in the framework of the MCM, including vrA^ and vrA states and also Pauli exclusion 



principle effects. Their conclusion was that data can be reproduced with a dipole form 
factor if A^rAf = 1.0 GeV and A^^a = 1.3 GeV. 

In this work we address inclusive meson production at high energy proton-proton 
and proton-nucleus collisions. We concentrate on low px {px < 0.3 GeV) and large xp 
[xp > 0.2) where non-perturbative effects are dominant and the meson cloud is most 
relevant. More specifically we address the xp distributions of vr"*" and K~^ produced in 
pp and pA collisions at piat — 100 GeV/c. We present simultaneously a fit of the meson 
inclusive spectra and the analysis of the GSR including the recently measured d{x) —u{x) 
distribution. Our purpose is to make a new test of the MCM and at the same time check 



the cut-off choices made in refs. 101 and 10 



2 Meson spectra in the MCM 

Assuming that in proton-proton collisions the proton behaves like a meson (M ) - baryon 
(B) state, the possible reaction mechanisms for meson production at large xp and small 
Pt are illustrated in Fig. 1. In Fig. la the baryon just "flies through", whereas the 
corresponding meson interacts inelastically producing a vr^ or a K^ in the flnal state. 
In Fig. lb the meson just "flies through", whereas the corresponding baryon interacts 
inelastically producing a 7?+ or a K~^ in the flnal state. In Fig. Ic the meson in the cloud is 
already a vr^ or a K~^ which escapes. We shall refer to the flrst two processes as "indirect 
production" and to the last one as "direct production" . The flrst two are calculated with 
convolution formulas whereas the last one is given basically by the meson momentum 
distribution in the cloud initial \MB > state. Direct production has been widely used 
in the context of the MCM and applied to study n, A"*""*" and 7r° production |0, 0, p!2| . 
Indirect meson production has been considered previously in a simplifled approach ||T3[| . 
In ref. a process analogous to Ic, in which the baryon escapes was used to determine 
various cut-off parameters. The authors found values around 1 GeV for all of them. 
However the agreement between theory and experiment was poor. This suggests that 



the no-rescattering assumption is not appropriate. In |T2[ the same type of process was 
considered but rescattering of the baryon with the target was included and absorptive 
corrections were calculated. This results in a reduction of about 40% of the cross section at 
almost all values of the baryon fractional momentum (z) except at z ~ 1 where absorptive 
corrections are not important. The indirect mechanism has not been applied to baryon 



production because the sub-process Bp -^ BX is not experimentally well known. On the 
other hand we can study inclusive meson production with la) since data on Mp -^ MX 
are available. 

As stated above, in the MCM the projectile proton is regarded as being a sum of 
virtual meson-baryon pairs and a proton-proton reaction can thus be viewed as reaction 
between the "constituent" mesons and baryons of one proton with the other proton. 

We shall decompose the proton in the following possible Fock states: 

\p> = Z[ Ipo > + bo7r° > + |r27r+ > +|A°7r+ > +|A+7r° > +|A++7r" > 

+ |Air+ > + |s°ir+ > + \TP*K^ > + |S+is:° > + |S+*i^° > ] (1) 

where \pq > is the bare proton. We consider only light mesons. In the state \Np > 
was also included but was found to be relevant only for production of large px particles. 
Since we will be restricted to low px production we shall neglect the havier mesons. 

The relative normalization of these states is in principle fixed once the cloud parame- 
ters are fixed. Z is a normalization constant. 

In the \MB > state the meson and baryon have fractional momentum i/m and ys with 
distributions called fM/MBiyu) and fB/MBiVB) respectively. Of course i/m + Vb = ^ and 
these distributions are related by: 

fM/MBiz) = fs/MBil - Z) (2) 

The "splitting function" fM/MB{y) represents the probability density to find a meson 
with momentum fraction y of the nucleon and is usually given by 

(3) 
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for baryons belonging to the decuplet. In the above equations t and uim are the four 
momentum square and the mass of the meson in the cloud state, t^ax is the maximum t 
given by 

i^ax = Mly - ^ (5) 



Mb (Mp) is the mass of the baryon B (p). Since the function fM/MB{y) has the interpre- 
tation of "flux" of mesons M inside the proton, the corresponding integral 

nM/MB = J dy fM/Msiv) (6) 

can be interpreted as the "number of mesons" in the proton or "number of mesons in the 
air" . In many works the magnitude of the multiplicities hm/mb has been considered a 
measure of the validity of MCM in the standard formulation with MB states. If these 
multiplicities turn out to be large ( ~ 1) then there is no justification for employing a 
one-meson truncation of the Fock expansion. The model has no longer convergence. This 
may happen for large cut-off values. 

The invariant cross section for production of positively charged mesons M"*" (vr^ or 
K~^) is given by: 



dp^ 71 dxpdp'j 



^ —3 = — .^ .^2 = ^M + $B + $i. (7) 



where 



and 



*- = ^X. jf-i-B{y) - ,^^^/^),^. My) (8) 

*- = St ^^^'^'^^'^ -ydixp/y)dpl^'^'y^ (^) 

$A/ and $B refer respectively to the indirect meson and baryon initiated reactions and xp 
and pt are respectively the fractional longitudinal momentum and transverse momentum 
of the outgoing meson. The sum is over all the cloud states in (|T|). 

$£) represents the direct process depicted in Fig. Ic and is given by: 

$D = E fM/MB{xp,pl) a^P(sx) Kabs (10) 

MB 

where fM/MB{xp,p'^) is given by (H) and (§) , not integrated over t and with the replace- 
ment t = — p|i/(l — Xp) — a;|'M^/(l — Xp). The quantity cr^^{sx) is the baryon-proton 
cross section at cms energy y/sx and Kabs is an absorption factor. This cross section is 
energy (sx) dependent but in the energy range considered here its variation is very small 
and therefore we take it as a constant cr^^ = 38 mb. The K factor was used in refs. 



12] , |T4| to account for rescattering of the escaping cloud element (in that case a baryon 



and here a meson) against the proton target. It varies in the range < Kabs < 1 and 



expresses the efficiency of the direct process. Of course this is a model dependent quan- 
tity. As stated above, the modest agreement of direct production calculations with data 



obtained in several works and the improvement obtained in [|l^, |Tj] (with the inclusion of 
absorption effects) strongly suggest that this factor is important and we shall keep it. 

Once the " splitting functions " (|^) and (^ are known we can calculate the part of the 
antiquark distribution in the proton coming from the pion cloud with the convolution: 

q{x)f= f'^fM/MB{y)rf{-) (11) 

Jx y ' y 

where q}{z) is the flavour / valence antiquark distribution in the pion. With the above 
formula we can compute the d and u distributions, their difference, D = d{x) — u{x), and 
calculate the Gottfried integral: 

Sg = 1-1 C\d{x)-u{x)]dx (12) 

O O JO 

3 Inputs for the calculation 

3.1 "Elementary cross sections" 

The invariant cross sections for the reactions Mp — > M^X and Bp — » M+X, appearing 
in the convolutions (P) and (^ are listed in Table I. Those in the first column can be 
taken directly from experimental data|T5|, |T^. The unmeasured vr^p -^ M^X process 



was approximated by the average between vr+p -^ M^X and vr^p -^ M^X . The same 
procedure was used for the process K^p —>■ M~^X. 

In making use of experimental cross sections we are treating the virtual mesons in the 
cloud as real particles. In our case this can be justified because the observed mesons both 
on the left hand side of (0) and on the right hand side of (^ and of {^ have very small 
tranverse momentum. Consequently the involved cloud particles must also have small 
tranverse momentum, being therefore not far from the mass shell. 

The reactions in the second and third columns of Table I are not measured. We expect 
them all to be of the same order of magnitude (aditive quark model approximation). We 
therefore approximate all these cross sections by an "average cross section " . In the case 
of TT"*" and K~^ production we shall assume respectively that: 

xpda^P^^^^ xpdaPP^^''^ xpda^P'^^''^ xpdaPP^^''^ 

r-^ - r-^ (13) 

Tidxpdpj^ -Kdxpdpj^ ' Tidxpdpjp Trdxpdpj^ 



This implies that all baryons are equally efficient as the proton to produce tt'^ or K^ in 
collisions with a target proton. The absorption factor appearing in (0) is chosen to be 
Kabs = 0.4 for pion production and Kabs = 0.8 for kaon production. These values are 
within the range of theoretical estimates presented in ^^ . 
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Table I: Subprocesses contributing to inclusive meson production (M+ = n^ or K^). 



3.2 Coupling constants 

The coupling constants are either measured, inferred from isospin symmetry or estimated 
with, for example, QCD sum rules. We will take them from other works P, 
keep them fixed. They are given in Table II. 
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Table II: Coupling constants. 



3.3 Cut-off parameters 

In the calculations we need the baryon-nieson-baryon form factors appearing in the split- 
ting functions. Following a phenomenological approach, we use the dipole form: 

FmbM=( ''T~ '"}' )' (14) 

In the above formula t^uBp is the form factor cut-off parameter. 

In [Q an extensive discussion concerning the appropriate value of the cut-off has been 
made, using exponential form factors. The conclusion of the authors was that, for ex- 
ponential form factors, A^^viv - 1000 ^eV, A^^va - §00 MeV and A^^y ~ 1200 MeV . 
They find out that allowing the tcNN and ttNA vertices to be different improves the 
quality of their fits. On the strange sector, on the other hand, they use a harder and 
universal cut-off. In ref. [0 the MCM was used to study baryon {n, A and A++) spectra 
in pp collisions and the resulting fits strongly suggest a universal cut off A*^ ~ 1000 MeV. 
The exponential form factor is, of course, not the only possible choice. One might use 
a monopole or dipole form as well. Differences between the various forms are not par- 
ticularly important and, besides, as it was pointed out by Kumano |T^ it is possible to 
"translate" one form factor with its cut-off to another form factor with a correspondingly 
different cut-off, the overall results being approximately equivalent. The approximate 
relation between the exponential (e), dipole (d) and monopole (m) cut-offs is given by 
A*" ~ 0.62A'^ ^ O.TSA*^. As pointed out in [§, in dealing with decuplet splitting functions 
monopole form factors lead to divergencies. With an exponential form factor one avoids 
this problem but the calculations can only be done numerically. Using the dipole form 
factor we can perform the integrations (§) and @ analytically. Because of this advantage 
we choose the dipole form. 

Once translated to the monopole form the cut-off values quoted above give values 
around 800 MeV which are still significantly smaller than old analyses performed with 
the Bonn potential [|l7l or the Nijmegen [|l9l potential which favour a harder cut-off 
(A™ ~ 1000 MeV) but are close to values obtained in more recent analyses performed by 
the same group pO|l . 

In view of all the works done so far on this subject, we may conclude that these cut-off 
parameters must be soft. The next question which is now under debate is: which cut-off is 
larger, Aj^nn or A^ttva ? As pointed out in the introduction, the E866 analysis [0 favours 



Attnn > A^TVA whereas Melnitchouk, Speth and Thomas [10| suggest that A^tva^ < A^rAfA- 



Inspired on these two works we shall test the following choices for the dipole cut-off 
parameters: 

(/) Kct = 0.96 GeV A,ec = 0.77 GeV 

{II) Kct = 0.87 GeV Kec = 1.0 GeV (15) 

Koct and A^ec are the cut-off parameters for all the octet and decuplet vertices respectively. 
As it will be seen, these choices are in reasonable agreement with experimentally measured 
meson xp spectra. 

4 Results and discussion 

Before presenting a comparison with experimental data we show in Fig. 2 the behaviour 
of the splitting functions fu/MB and fs/MB with respect to cut-off variations. Solid and 
dashed lines represent respectively the octet eq.(^ and decuplet eq.(|) MB states. The 
two upper (lower) boxes show functions calculated with soft (hard) cut-off parameters. 
The upper and lower boxes on the left (right) show functions with octet cut-off parameters 
harder (softer) than the decuplet ones. The values, chosen for the purpose of illustration, 
are: Aoct = 1.0 GeV and Adec = 0.80 GeV (Fig. 2a); A^^ = 0.80 GeV and Adec = 1.0 GeV 
(Fig. 2b); Aoct = 2.0 GeV and Adec = 1.6 GeV (Fig. 2c); Aoct = 1.6 GeV and Adec = 2.0 
GeV (Fig. 2d). We consider the momentum (y) distributions of vr^ in the states vr+n 
and 7r"*"A° and also the baryon momentum [ys = 1 — y) distribution in these states. 
Because of eq.(^) all these figures are symmetric and the meson momentum distributions 
(always on the left part of each box) are the "mirror" pictures of the corresponding baryon 
momentum distributions (always on the right part of each box). The comparison between 
the upper pannels with the lower ones shows that at soft cut-off's there is a more distinct 
separation between baryons and mesons, the former peaking at higher values of y and 
the latter at smaller y's. At harder cut-off's all the curves tend to become identical and 
centered at y = 0.5. The decuplet functions lie always bellow the octet functions, except 
if Adec ^ Aoct, which is unrealistic. In Fig. 2d we can see that f-K+fn+Aoiy) lies above 



/j^+/^+„ for y > 0.6. This feature was found in ||10[ to be very important in the study of 
sea antiquark asymmetry. 

Using, as input, the splitting functions shown in Figs. 2a, 2b, 2c and 2d in eq. (|^ 
for M"*" = 7r+ we obtain the pion spectra shown in Figs. 3a, 3b, 3c and 3d respectively. 

9 



They are decomposed in meson initiated (M), eq. (^ (dot-dashed hnes), baryon initiated 
(-B), eq.(P) (dashed hnes), direct (D), eq.(0) (dotted hnes), and the total sum (T). 

We see that the direct process dominates pion production up to xp — 0.7. In the 
two lower figures (harder cut-off parameters) they are so important for all xp that one 
can neglect the M and B contributions. For soft cut-off choices the B component be- 
comes important too and even dominant at xp > 0.8 (shaded area in Figs. 3a and 3b). 
Comparing the two shaded areas, we observe that with AjrAfAr < A^jva (on the right) 
the B component starts earlier to dominate the pion xp spectrum. This component is 
the harder and flatter one and is essential to obtain a good description of data at large 
xp. The direct contribution alone would fah too fast and underestimate data. Neglecting 
the direct contribution leads to a good flt only with the M and B components but at 
the expense of unrealist large (dipole) cut off choices: Aod = 1-13 GeV and Adec = 1-8 
GeV. Only with these large numbers we recover the normalization of the experimental 
spectrum. Neglecting also the B component leads to a good flt with the M component 
alone and a cut-off Aod = 2.1 GeV. In short, it is clear that we need all three components 
to describe data, the flrst one (D) to get the proper normalization and the last two {B 
and M) to get the correct shape of the meson spectra. In this way we get also reasonable 
pion (and also kaon) multiplicities. For both cases (I and II) we obtain n^r — 0.17, which 
is small enough to justify the use of the one- meson truncation of the Fock decomposition. 

In Fig. 4a and 4b we present our pion and kaon spectra compared with experimental 



data, circles from ref. [|I5| and triangles from ref. [|I6|. In these flgures we include all 
three components and plot eq.(|^ with the cut-off choices I (solid line) and II (dashed 
line). As it can be seen, the overall agreement with data is good over a large range on xp. 
We have checked that signiflcant deviations from these choices lead to large disagreement 
between our spectra and data. As an example we show in Fig. 4 with dotted lines the 
spectra corresponding to the the choice Aoct = 1.0 GeV and Adec = 1.3 GeV discussed in 
[ p!0| . As it can be seen they are already somewhat far from data points, especially in the 
case of the kaon spectrum. From this flgure we may conclude that the MCM can give 
an accurate description of data with soft (A < 1.0 GeV) cut-off parameters. Choosing 
Aoct > Adec or Aoct < Adec docs not make much difference as long as both cut-off's differ 
only by ~ 200 MeV. 

In Figs. 5a and 5b we show the difference d — u and ratio d/u respectively. The quark 
distribution functions were calculated with the convolution formula (|TT|). The points 
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represent the E866 data. Solid and dashed hnes represent the parameter choices I and II 
respectively. As expected our solid lines are very close to the ones presented by the E866 
analysis in 0| . Using a larger cut-off for the ttNA vertices leads to the dashed curves. As 
pointed out in [|T0| this choice leads to a reduction in the antiquark asymmetry, leaving 
room for additional sources of asymmetry, such as the Pauli exclusion principle. 

The value of the Gottfried integral given by (|l^) is Sq = 0.261 (case I) and Sc = 0.314 
(case II). The value reported by the NMC collaboration is 0.235 ± 0.026. Taken together. 
Figs. 4 and 5 and the integral 5*^, it seems that choice I is better everywhere but fails 
badly in the d/u ratio. Choice II, on the other hand, is never far from data points but 
gives a too large value for So- We thus conclude that a better global description of data 
can be obtained with A^d moderately smaller than A^ec and that some additional source 
of asymmetry, other than the meson cloud, is required. 

Another interesting test for the MCM is inclusive meson production in proton nucleus 
coUisions. In principle we can use the same formulas eqs. (|^, (^, (y) and (jT^), replacing 
the cross sections cr^^(-B)+p^Af x -^^ ^m{b)+a^m x ^^^ keeping exactly the same splitting 
functions. However most of the M{B) + A —>■ M^X cross sections are not measured. As 



a first test we shall assume, following the experimental analysis described in refs. [|^, [16 
that: 

j3 pA^M+X j3 pp^M+X 

E—-—- ~ const A" E—-—- (16) 

dp-^ dp-^ 

where the constants are chosen to reproduce the normalization of data points. 

Usually a is a very slowly varying function of xp- We take it constant {a = 0.72). 

The main goal of this exploratory study is only to test the shape of the spectra obtained 

with the MCM. In the above expression the pp cross sections are exactly the same as 

before. In Fig. 6 we show the spectra obtained with (|16D for pions (6a and 6b) and for 

kaons (6c and 6d). The cut-off choice I (II) is depicted with solid (dashed) line. Since the 

\MB > state will interact inside the nucleus it may be affected by medium effects. The 

simplest way to incorporate these effects is by changing the baryon masses according to 

Mb —>■ M^ ~ 0.8Mb. With this modification the meson splitting functions will increase 

in normalization and will peak at larger values of yM- Because of this the final meson 

spectra will also peak at larger values of xp. The inclusion of medium effects (curves in 

dotted lines in Fig. 6) improves the agreement between the MCM and data. 
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5 Conclusions 

This work was motivated by the recent revival of the meson cloud picture of the nucleon. 
We performed a comparison between inclusive meson spectra predicted by the MCM and 
experimental data. Instead of using only the "direct process" , as done in refs. [0, 0] we 
have included the "indirect process", in which the final cross section is a convolution of 
the MCM "splitting functions" with the "elementary" meson-nucleon and baryon-nucleon 
cross sections, in an analogous way as done in the QCD parton model calculations. In this 
way we can test the universal splitting functions and determine the cut-off parameters. 
This determination is done by simultaneouly analysing mesonic spectra, difference and 
ratio of antiquark distributions and the Gottfried sum rule. 

The first conclusion of this work is that the MCM describes reasonably well the pro- 
duction of pions and kaons in pp high energy collisions. Even the intermediate xp region 
(0.2 < xp < 0.6) is reproduced. 

Concerning the relative strength of the nucleon and delta vertices, we find that the 
MCM can give a good description of data with soft (A < 1.0 GeV) cut-off parameters. 
Choosing Aoct > ^dec or Aoct < A^ec does not make much difference as long as both cut-offs 
differ only by ^ 200 MeV. The latter choice seems to be more appropriate, especially if 
one introduces additional sources of sea asymmetry. 

Acknowledgements: This work has been supported by CNPq and FAPESP under 
contract number 98/2249-4. We are indebted to Y. Hama, P.M. Steffens, A.W. Thomas 
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Figure Captions 

Fig. 1 pp or pA collision in which the projectile is in a \MB > state, a) "indirect" K~^ 
or 7r+ production: the meson M from the cloud undergoes the reaction M p{A) — > 
7i^{K^)X. b) "indirect" fC^ or tt^ production: the baryon from the cloud undergoes 
the reaction B p{A) -^ 7r+(K+)X. c) "direct" K^ or 7r+ production: the cloud 
meson vr^ (-^^) escapes from the cloud as a spectator. 

Fig. 2 Illustration of the octet (decuplet) splitting functions /a//j\/b and fs/MB in solid 
(dashed) lines for the following cut-off choices: a) Aoct = 1.0 GeV and K^ec = 0.80 
GeV; b) Aoct = 0.80 GeV and Adec = l-O GeV; c) Aoct = 2.0 Ge\/ and A^ec = 1.6 
Gel^; d) A^,* = 1.6 GeV and A^ec = 2.0 Gel/. 

Fig. 3 Total (T) inclusive pion spectra in pA collisions calculated with eq. (0) and the 
direct (-D), meson initiated (M) and baryon initiated (B), given respectively by 
eqs.([TO|), (H) and @. a), b), c) and d) correspond to the cut-off choices made in 
Fig. 2. 



Fig. 4 a) Inclusive pion spectra calculated with (|^. Data are from |]T5|, |T6|. Solid, dashed 



and dotted lines correspond to the cut-off choices I, II and A,rAf = 1.0 GeV, A^^a 
1.3 GeV respectively; b) the same as a) for kaon spectra. 



Fig. 5 a) d ~ u calculated with eq. ([TT|) with cut-off choices I (solid line) and II (dashed 
line) compared with E866 data; b) the same as a) for the ratio d/u. 

Fig. 6 Inclusive meson spectra in proton nucleus collisions calculated with ([T6|) using cut- 
off choices I (solid lines) and II (dashed lines). Dotted lines show the curves obtained 
taking medium effects into account. Data are from |l^, |T^. a) pion spectrum with 
choice I; b) pion spectrum with choice II; c) same as a) for kaons; d) same as b) for 
kaons. 
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